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Abstract. We present theoretical predictions of the rotational line emission of deuterated water in low-mass 
protostar collapsing envelopes. The model accounts for the density and temperature structure of the envelope, 
according the inside-out collapse framework. The deuterated water abundance profile is approximated by a step 
function, with a low value in the cold outer envelope and a higher value in the inner envelope where the grain 
mantles evaporate. The two abundances are the two main parameters of the modeling, along with the temperature 
at which the mantles evaporate. We report line flux predictions for a 30 and 5 Lq source luminosity respectively. 
We show that ground based observations are capable to constrain the three parameters of the model in the 
case of bright low-mass protostars (L>10 Lq), and that no space based observations, like for example HSO 
observations, are required in this case. On the contrary, we show that the study of low-luminosity sources (L<10 
L©), assuming the same HDO abundance profile, requires too much integration time to be carried out either with 
available ground-based telescopes or with the HIFI instrument on board HSO. For these sources, only the large 
interferometer ALMA will allow to constrain the HDO abundance. 
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1. Introduction 

Water is a key molecule in the study of star formation for 
four main reasons. Since it has a large dipole, it has large 
spontaneous emission coefficients and high critical densi- 
ties, which makes it a powerful diagnostic tool. Secondly, 
water is important in star f orming regions for the rol e 
it plays in cooling the gas llNeufeld fc Kaufman! 1199.1) 
and, sometimes, also in heating it (e.g. ICeccarelli et alJ 
Il996l hereinafter CTH96). Third, it is also an impor- 
tant molecule from the chemical point of view, because 
in many circumstances it is the most abundant oxygen- 
bearing molecule, and, hence, it influences the abundance 
of all the more complex "trace" molecules; in other words, 
it infl uences the overall chemical composition of the gas 
(e.g. iRodgers fc Charnlevl I20031. Finally, water is cer- 
tainly the most abundan t component of grain mantles too 
(e.g. iGibb et alJ 12004^1 . For all these reasons, it is ex- 
tremely important to know its abundance distribution in 
star forming regions. Unfortunately, the interstellar ro- 
tational water lines are absorbed by the water vapour 
of the Earth atmosphere. Therefore, the observation of 
the water vapour in astrophysical sources requires out-of- 
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atmosphere instruments, like those on the past satellites 
ISO, SWAS and ODIN, or the future satellite HSO, to be 
launched in 2007. 

Fortunately enough, some rotational transitions of the 
deuterated form of water, HDO, are observable from 
ground based telescopes, and can be used, under some 
circumstances, to probe the water content in astrophysi- 
cal sources. Furthermore, HDO is an interesting molecule 
on its own, because it gives important information on the 
water formation route on the one hand, and, on the other 
hand, on the physical conditions of the observed region, 
either at present or during the past. If formed on the grain 
surfaces during a previous phase, for example, and then 
released into the gas phase because of the mantle evapora- 
tion due to the heating from a new born star, the measured 
HDO abundance would bring the hallomark of the past 
history of water formation. This is likely the case of what 
happens in solar type protostars. After a pre-collapse cold 
phase where the grain mantles are likely formed, either 
because of accretion of molecules onto the grain surfaces 
or because of the synthesis of the molecules on the grain 
surfaces, the molecules are released in the gas phase in the 
regions heated by the protostar radiation. Two classes of 
molecules show evidence of this process: 
i) Hydrog enated molecules like formaldehyde and 
methanol (|Ceccarelli et alJ l2000bl ISchoier et ail I2002L 
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l2004t iMaret et all 12004 |2005[) and complex organic 
molecules, observed in the so called "hot cores" 
llCazaux et alJl2003t iBottinelli et~aill2004alrJ: iKuan et all 

2004J. These molecules are not observed in the cold gas 
with the large abundances measured in the very inner re- 
gions of the envelopes surrounding solar type protostars. 
Indeed, they are direct and indirect products of the grain 
mantle evaporation, so that they are, in a way or another, 
linked to the pre-collapse phase. 

ii) Deuterated molecules observed with extremely large 
abundances, which can only be formed during the 
cold, de nse and CO depleted phase of the pre - 
collapse l|Roberts fc Millarl l2000t iRoberts et all 120031) 
Nota ble exampl es are formaldehyde jCeccarelh' et alJ 
Il998l) . methanol jParise et alJl2002l l2004|h hydrogenated 
sulfide JVastel et al]l2003h and ammonia l|Lis et alJl2002l 
Ivan der Tak et alJl2002V 

HDO belongs to the second class of molecules witnesses of 
the past pre-collapse phase. And, it is, apparently, a very 
different beast from that zoo. Indeed, the above mentioned 
molecules show extremely large deuteration ratios in low 
mass star forming regions: singly deuterated molecules 
are more than 10% of their H-counterparts, whereas dou- 
bly and triply deuterated molecules can respectively be 
as high as 10% and ~ 1%. On the contrary, HDO is 
only a relatively sm all fraction of H2O, about 3% at best 
llParise et al.l l2005h . and no D2O has been detected so 



far, at relatively low ratios to our knowledge (Cernicharo, 
private communication). The above relative small values 
agree indeed also with the non-detec t ion of solid HDO i n 
protostars (e.g. iDartois et aT]|2003t IParise et alJl2003|) . 
Thus, there is evidence that the water deuteration fol- 
lows a different route with respect to the other molecules. 
In addition, the study of the available observations in 
IRAS16293-2422 suggests that the deuteration ratio is not 
the same in the inner and outer envelope: 3 % and <0.2% 
in th e inner and outer envelope respectively l|Parise et alJ 
2005). Thus, in order to better understand this difference, 
one needs to be able to disentangle where the different 
emission comes from. The present study aims at giving a 
theoretical tool for the estimate of the HDO abundance 
in protostellar envelopes, in the inner and outer regions, 
and, possibly, the evaporation temperature, and, conse- 
quently, at guiding the best observations to obtain those 
informations. Specifically, given that the satellite HSO will 
be launched in a relatively near future (2007), a particular 
care will be devoted to discuss whether observations with 
instruments on this satellite will be required for deriving 
the HDO abundance and evaporation temperature in low 
mass protostellar envelopes, against the observations that 
can be obtained with ground based telescopes. 

The article is organized as follows: we describe the 
adopted model and code in §2, we present the HDO 
molecule as well as the instruments available to observe its 
rotational spectrum in §3. We describe the general charac- 
teristic of the predicted spectrum in §4 and give the results 
in §5 for a relatively high (30 L ) and low (5 L©) lumi- 
nosity case respectively. Finally, §6 contains a discussion 



about the best lines to observe to derive the HDO abun- 
dance in the inner and outer envelope, as well as the HDO 
ices evaporation temperature. We also discuss the interest 
of carrying space based versus ground based HDO obser- 
vations to constrain these parameters. 

2. The model description 

The model used to compute the HDO emission in collaps- 
ing envelopes has been adapted from the original model 
developed by CHT96 to predict the OI, CO and H 2 
line emission, and s uccessively modified to compute the 
H2CO line emission (jCeccarelli et al.ll2003fl . Here we de- 
scribe very briefly the basic assumptions of the model. 

The envelope density structure and dynamics ar e as- 
sume d to follow the "inside-out" collapse picture ijShvJ 
Il977h . for a spherical initial isothermal sphere undergo- 
ing collapse. In the outer envelope, not affected by the 
collapse yet, the molecular hydrogen number density dis- 
tribution nn 2 (f) is given by: 



riH 2 (r) = 



2-KfmiHG 



(1) 
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where a is the sound speed, uih is the hydrogen mass, /1 
is the mean molecular mass in amu units, equal to 2.8, 
fiooAU is the distance from the center in 100 AU units, 
and G is the gravitational constant. In the inner collapsing 
regions the density is described by the free-fall solution: 
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The free fall velocity is given: 
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where M is the mass accretion rate, related to the sound 
speed by 

. 3 

M = 0.975^- = 10~ 5 ' 

G V 0.35km s 



(035=7*) Msyr " (4 » 



M_5 is M in units of 10 -5 M Q yr -1 , and M*i is the mass of 
the central object M* in units of 1 M Q . The gravitational 
energy of the collapsed mass is released radiatively as ma- 
terial falls onto the core radius i?* , so that the luminosity 
of the central object is: 



GM*M ( M± 1 • r , T 
where R12 =R*/10 12 cm. 



(5) 
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The spherical symmetry is conserved through the col- 
lapse in this model. In this sense, the model gives accu- 
rate results only for that part of the envelope in which the 
spherical symmetry is a good approximati on, i.e. probably 
for rad ii larger than a few tens of AUs l|Ceccarelli et alJ 
20003- At smaller scales large deviations from the spher- 
ical symmetry are expected because of the presence of a 
circumstellar disk. 

The gas temperature is computed self-consistently by 
equating the cooling and heating mechanisms at each 
point of the envelope. Details of these calculations can be 
found in CHT96. What is important to notice here is that 
the cooling depends on the abundance of the atomic oxy- 
gen, CO and H 2 0, namely the main gas coolants across 
the envelope. Hence their abundances are hidden parame- 
ters of the model. The resulting gas temperature tracks the 
dust temperature very closely across the envelope, except 
where the icy mantles evaporate. Because of the injection 
of large quantities of water in the gas phase and the con- 
sequent enhancement of the gas cooling, the gas and dust 
(thermally) decouple in a relatively small reg ion after the 
icy m antles evaporation (see also for example iMaret et alJ 
l2002jL 

ICeccarelli etaH l|2000al) and IMaret et all l|2002|) 
demonstrated that the gas temperature differs by no 
more than 10% from the dust temperature, and this oc- 
curs in a relatively small region just next to where the 
ices sublimate. There is some uncertainty in the water 
and oxygen abundances in protostars, for derivations of 
the H 2 abundance has been obtained in two sources 
only, IRAS 16293 and N GC1333-IRAS4 l|Ceccarelli et all 
2000a HMaret et al.l2002ft . However, these uncertainties do 
not affect much the gas temperature derivation, for in the 
outer envelope the gas is cooled mostly by CO 1 , and the 
inner is cooled by the sublimated water. At the same time, 
in the inner region the heating of the gas is dominated by 
the collisional de-excitation of the water molecules, pho- 
topumped by the photons emitted by the innermost warm 
dust (CHT96), so that, in the end, the coupling between 
the dust and gas is always ensured. As a result, the HDO 
computed line intensity are rather insensitive to these un- 
certainties. 

The original CHT96 model is time-dependent, but the 
authors demonstrated that the physical structure and wa- 
ter chemistry can be computed at each time independently 
on the previous history, once the luminosity and the mass 
of the protostar are known. The reason for that is the fol- 
lowing. In the outer cold envelope the abundance of the 
water is practically constant with time, whereas in the in- 
ner warm envelope the water abundance is set by the ice 
evaporation, except in the very inner regions where the 
temperature exceeds 250 K. Since the sizes of the evapo- 
ration region are set by the temperature and the density 
profile, for practical applications to real sources one needs 
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to input them only. In addition to that, the velocity field 
is also needed to compute the line intensity, and this can 
be derived as well by observations. 

For the standard case described in section 4, we consid- 
ered the case of IRAS 16293, whose velocity, density and 
temperature profiles have been derived by several authors 
from actual observations. The co mputations show n in Fig 
!2l7l adopted the profile derived in lCeccarelli et akl l|2000a[) . 
We also studied the case of a low luminosity source, and, at 
this scope, we used the velocity, density and temperature 
profile s derived for the source L1448mm bv lMaret et alJ 

Finally, we approximate the HDO abundance profile 
with a step function. In the outer envelope the HDO abun- 
dance, x co id, is relatively low and similar to that observed 
in molecular clouds. In the innermost regions of the enve- 
lope, where the dust reaches the evaporation temperature 
Tevap, the HDO abundance, x warm , jumps to larger values. 
These last three parameters — x co \d, £ W arm and T cvap - 
are the three free parameters of the present modeling, and 
they will be varied to study their influence on the HDO 
line emission. 

Low-mass protostars exhibit energetic molecular out- 
flows, where the gas can be very hot (up to 2000 K, as 
testified by the H 2 rovibrational emission observed in sev- 
eral outflows). This gas could in principle contribute to 
the HDO line emission, so that it is important to evaluate 
whether the outflow emission can be comparable or even 
larger than the HDO emission from the envelope. From a 
theoretical point of view, since molecular deuteration de- 
creases exponentially with increasing gas temperature, the 
HDO abundance is expected to be v ery low in the warm 
gas of molecular outflows (see also iBergin et al1ll999l) . 
To support these theoretical expectations, observations of 
outflows on several low-mass protostars have revealed no 
HDO emission, at a level at least 10 times fainte r than 
towards the protostar (TRAS1 6293. IParise et all 120051: 
NGC1333-IRAS2 and IRAS4A, ICaux fc a.lJlin nrenh . We 
therefore decided to consider in the present study only the 
envelope contribution to the HDO emission. 

All the line fluxes are computed for a source 
at a distance of 160 pc. The model computes the 
level population of the first 34 levels of the HDO. 
The molecular data are from the JPL catalogue 
(http://specJpl.nasa.aov/home.html), and the collisional 
coefficients are from iGreenl ljl989j) . for He-HDO collisions 
between 50 and 500 K, scaled for collision with H 2 . Note 
that all abundances are reported here with respect to H 2 . 

3. The HDO rotational line spectrum 

3.1. The energy levels 

HDO is a planar assymetric top molecule. Its dipole lies 
along its a and b axies . with dipole mome nts \i a = 0.657 
and fj, b = 1.732 Debye ijcioueh et alJll973|) . The allowed 
transitions are thus a-type and 6-type transitions : 
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and AK, = ±1 
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or Receiver), will observe in the 500-2100 GHz band 

. T , . . i a ii and GREAT (German REceiver for Astronomy at 

= ±l and AK C = ±l . freque l ncies) in ^ l>w THz 

The first energy levels for HDO are presented in Fig.^ c) The HSO satellite 

The HIFI (Heterodyne Instrument for the Far- 
Infrared) instrument on board Herschel Space 
Observatory (launch planned in 2007) will allow ob- 
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4. Model results 

4, 4 ^=^ = 3 2 ,i 

50 - 40 4 3 2 ,2 We first report the computed HDO line spectrum of a 

test case, to illustrate the general characteristics of the 
predicted spectrum. Then, in the following paragraphs we 
discuss thoroughly how the HDO line spectrum varies with 
the three free parameters of the model, namely the HDO 
abundance in the outer cold envelope x C oid, and in the 
inner warm envelope x waTm , and the evaporation temper- 
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Fig. 1. First energy levels for the HDO molecule. 

The list of the transitions between the first 34 levels is 
reported in Table IA~T1 

3.2. The observing instruments 

We detail in this section the instruments that operate (or 
will operate in the near-future) in the frequency bands 
interesting for the HDO observations. 

a) Ground telescopes 

Several ground based facilities allow to observe the 
HDO rotational spectrum transitions. 
The IRAM (Institute for Millimetre Radioastronomy) 
30 m telescope (Pico Veleta, Spain) operates between 80 
and 281 GHz. 

The 15 m JCMT dish ( James Clerk Maxwell Telescope, 
Mauna Kea, Hawaii) observes in the 230, 345, 450, and 
650 GHz frequency bands. The 800 GHz instrument is no 
longer available, but this frequency range will be observ- 
able with APEX. 

The 10 m CSO dish (Caltech Submillimeter Observatory, 
Mauna Kea, Hawaii) operates roughly in the same fre- 
quency bands as JCMT. 

The 15 m APEX dish (Atacama Pathfinder Experiment, 
Chajnantor, Chile) is being built on the ALMA site. It 
will allow observations in the 230, 345, 650 and 850 GHz 
bands. 

b) Airborne instruments: SOFIA 

The 2.5 m telescope of the Stratospheric Observatory 
for Infrared Astronomy, on board a Boeing 747, will 
be equipped with several interesting instruments for 
HDO observation. In particular, CASIMIR (CAltech 
Submillimeter Interstellar Medium Investigations 



Table 1. Values of the parameters adopted for the test 
case. The upper part reports the values of the three free 
parameters of the model, whereas the lower part reports 
the value of the hidden parameters. 



Parameter 


Value 


^warm 


1 x 10" 7 


^cold 


1.5 x 1(T 10 


T C vap 


100 K 


M* 


0.8 M 


M 


3 x 10~ 5 Meyr -1 


x(CO) 


1 x 10 -4 


x(O) 


2.5 x 10~ 4 


X C old(H 2 0) 


5 x 10~ 7 


Xwarm(H2 0) 


3 x 10~ 6 



The parameters adopted in the test case are reported 
in Tabled The values of the hidd en parameters are the 
same of the test case shown in ICeccarelli etall l|2003|) for 
discussing the H2CO line emission, and have been derived 
from several observations of IRAS16293-2422. The x C oid, 
£warm and T evap parameters have also been derived by 
the an alysis of the HDO o bservations towards IRAS16293- 
2422 (|Parise et al.ll2005|) . In the following all fluxes are 
given in ergs -1 cm -2 and are the integrated emission over 
the entire envelope for a source at 160 parsecs. Although 
the conversion in Kkms -1 would have been more conve- 
nient for the observers, the signal in Kkms -1 depends on 
the beam of the telescope used for the observation and 
how it matches the predicted extent of the emission. We 
will give the Kkms -1 signal for a few lines, when consid- 
ering the application to observations with ground versus 
space based instruments (§5). 

Table lA~T1 reDorts the predicted fluxes of the HDO lines 
for the test case, and Fig. |21 shows the line fluxes as a 
function of the upper energy of the transition. 
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Fig. 2. HDO line fluxes predicted for the test case 
IRAS16293— 2422 versus the upper level energy of the 
transition. The upper panel presents the transitions ob- 
servable from the ground, the middle panel the transi- 
tions observable with HSO/HIFI and the lower panel the 
transitions observable with SOFIA. The frequencies of the 
transitions that are going to be studied in more detail are 
indicated in GHz. 
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Fig. 3. Contribution to the line intensity as function of 
the distance from the central source, for seven HDO lines 
observable from the ground (upper panel) and six from 
space based telescopes (lower panel). 

to 150 cm -1 (~200K). This figure illustrates that the rela- 
tive contribution to the line intensity from the inner region 
(r < 150 AU) increases with increasing upper level energy 
of the transition. While all lines intensities are dominated 
by the emission in the inner region (where the abundance 
is about 3 orders of magnitude higher), the lines with low 
upper level energy exhibit a relatively high contribution 
from the cold outer envelope. 



The transitions have been labeled according to the 
available observing facility. It should be noted that the 
lowest lying levels are observable from the ground and 
that space borne observations become competitive for up- 
per energies greater than 50 cm -1 . 

Fig-Elshows the contribution to the emission from dif- 
ferent shells of the envelope for several transitions with dif- 
ferent upper energy levels, observable with ground-based 
telescopes (upper panel) or with Herschel-HIFI (lower 
panel). Note that the 893 GHz line has been included in 
the upper panel as it will be observable with the APEX 
facility. Some of the lines have been observed towards 
IRAS16293-2422 using the JCMT a nd the IRAM 30 m 
telescopes (80, 225, 241 and 464 GHz. IParise et alJl2005|) . 
The other ones have been chosen because they are poten- 
tially bright lines (Fig.0), spanning upper energy levels up 



4.1. Varying the parameters 

In this paragraph, we show how the HDO line spectrum 
changes varying the three free parameters of the model, 
namely the HDO abundance in the outer and inner parts 
of the envelope, x co \a and a; warm respectively, and the man- 
tle evaporation temperature T evap . 

^coid o,nd x warm 

Fig.0]shows the line intensities of eight HDO lines (in- 
cluding two lines observable with Herschel-HIFI) for a; co id 
varying from 10~ 7 to 10~ 12 and x warm varying from 10 -5 
to 10~ 10 . As expected, the intensities of the lower level 
transitions are mostly sensitive to the emission from the 
cold part of the envelope, as long as the outer abundance 
of HDO is not too low compared to the inner abundance. 
For instance the 464 GHz line is sensitive to a; co id even for 
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Fig. 4. Line intensity as a function of the HDO abun- 
dance in the outer and inner regions of the envelope, a; co id 
and x waxm respectively, for the eight HDO transitions at 
464, 893, 509, 80, 241, 225, 919 and 995 GHz respectively. 
The first six lines are observable from the ground, whereas 
the last two ones will be observable with the Herschel tele- 
scope. Line fluxes are in ergs^ 1 cm -2 . In these calculations 
the mantle evaporation temperature is 100 K. 

inner abundances up to two orders of magnitude larger. 
On the contrary, the higher the upper energy, the more 
the transition becomes dependent on the inner parts of 
the envelope. 

The figures also show that the low lying lines are much 
less sensitive in general to the variation of x W arm (a factor 
50 at most varying x warm by four orders of magnitude) 
with respect to the high lying lines. The latter vary almost 




log(X worm ) log(X worm ) 

Fig. 5. Predicted flux for the eight transitions depending 
on the assumed evaporation and inner HDO abundance 
for an outer abundance a; co id fixed to 2xl0 -10 . 



linearly with x W arm, so they are clearly better suited to 
constrain this parameter than the low lying lines. On the 
other hand, x co id is better constrained by the low lying 
lines, and particularly by the 464 GHz one, which shows 
the largest variation, nearly 2 orders of magnitude for a 
change of x co id of two orders of magnitude. 

^/ -^cvap 

Fig. [S] and [5] present the predicted fluxes depending 
on the evaporation temperature and respectively the in- 
ner and the outer abundances. The higher the evaporation 
temperature, the higher the inner abundance must be to 
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Fig. 6. Predicted flux for the eight transitions depending 
on the assumed evaporation and outer HDO abundance 
for an inner abundance £ warm fixed to 10~ 7 . 



reproduce the same flux. Indeed, the higher the evapora- 
tion temperature, the less extended the inner region. 

The 225 GHz line is mostly independant of the evap- 
oration temperature in the studied range (Fig O ■ Indeed 
this line has the highest upper energy (~ 170 K) and is 
likely to be mostly excited in the innermost region of the 
envelope. Thus it does not depend on the evaporation tem- 
perature, as long as this temperature is much lower than 
170 K, nor on a; co id- This line is thus only sensitive to the 
inner abundance. 

Another way to understand the dependance of the 
emission versus the evaporation temperature is to plot the 
regions spanned when varying the evaporation tempera- 
ture for a fixed emission flux in a ^coid-^warm diagram 
(Fig.0. 

As expected, we can notice that for low inner HDO 
abundances, the effect of the evaporation temperature is 
negligible, but curves begin to separate for high x war m- 
The uncertainty on the evaporation temperature mostly 
translates into an uncertainty in x warm . A remarquable 




log(X„ rm ) 



Fig. 7. Evaporation temperature dependence in a x warm - 
a^coid diagram. For each transition, the three sheaves cor- 
respond to the contour for three different values of the 
line flux. The different curves for each sheaf correspond 
to different values of the evaporation temperature. A nar- 
row sheaf implies that the uncertainty on the evaporation 
temperature value will not affect much the estimation of 

•£warm and £cold- 



point here is again that the 225 GHz emission is indepen- 
dent of the temperature. 

Fig. [5] shows that the line fluxes depend much more 
on x warm than on T eV ap- Actually, based on Fig. and El 
it seems very difficult to find a way to constrain T OV ap- 
However, Fig. 7 sheds some lights, and gives a sense of 
how the two parameters can be likely constrained. The 225 
GHz line flux in practice can be used to estimate x warm 
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(Fig. El and 0), irrespectively of T evap , whereas one can 
use, for example, the 464, the 80 and/or the 509 GHz 
line fluxes to estimate T cvap , knowing Xwarm, from Fig. [7] 
Evidently, it is not possible to easily disentangle the three 
parameters, but the above analysis suggests that a full 
modeling can possibly succeed. 



5. Applications 

5.1. IRAS 16293-2422 

In this section we try to give some practical recipes to 
derive, in first approximation, the values of the parame- 
ters of the model from actual observations. However, we 
emphasize that the described method will just give some 
approximative estimate of the parameter values, and that 
a full modeling of the source is necessary to derive more 
precise estimates. 

As a practical example , we will use the data obtained 
towar ds IRAS16293-2422 jstark et all 12004 IParise et all 
12005]) and we will derive the approximate values of x co id 
and x warm . We will then compare these results with those 
obtained with full mo deling the HDO em ission from this 
source, as described in lParise et alJ l)2005|) . 

First, to facilitate the comparison, Table El reports 
the conversion factor of the line intensity in erg s _1 
cm -2 into the velocity integrated main beam temper- 
ature T mb Av observed at the IRAM and JCMT tele- 
scopes respectivel y, for the five line s observed towards 
IRAS16293-2422 ijParise et alJl2005|) . Note that this is 
simply a conversion table, which does not take into account 
the convolution of the predicted emission with the telescope 
beam. This is strictly valid for a point-like source with re- 
spect to the telescope beam sizes. The reader should be 
aware that the convolution with the beam can introduce 
large factors of difference, if the emission originates in a 
region larger than that encompassed by the used beam, 
as it may be the case for low lying lines. In those cases, in 
order to use the plots of this article, the observer should 
take care to integrate the emission over the entire emitting 
re gion. 

IParise et alJ l)2005|) performed a full analysis of the 
HDO emission towards IRAS16293-2422. We describe 
however how the tools we have given here allow to ap- 
proximately derive the envelope parameters. According to 
Figure [7| the emission of the 225 GHz line, as well as the 
upper limit on the x C oid derived from the 464 GHz line, 
do not depend on the evaporation temperature. We can 
then use the Figure 01 (computed for an evaporation tem- 
perature of 100 K ) to derive an uppe r limit on a: co id from 
the 464 GHz line. IParise et al.l l|2005|) measured a line in- 
tensity of 5.5Kkms _1 . Using the conversion factor given 
in Table |U this corresponds to 2.2xl0 -15 ergs _1 cm -2 . 
We deduce from Figure 0] that x co id is lower than a few 
10~ 10 to 10" 9 . The 225 GHz line emission then allows 
to constrain the value of x waim (independantly of the 
evaporation temp erature). Ind e ed, co nverting the inten- 
sity observed bv IParise et all l|2005l) . we find 7x10" 17 



ergs _1 cm -2 , which with the condition x co id lower than 
10~ 9 , constrains the abundance of HDO in the warm en- 
velope to a value around 10 . The full modelling indeed 
leads to a; waTrn =(lzfc0.3)x 10~ 7 , a nd to a 3cr upper limit on 
x cold of 10~ 9 l|Parise et al.ll2005[) . 

Table compares the predicted fluxes, derived from 
the model by both using the simple conversion flux for a 
point-like source and an accurate convolution with the 
telescope (IRAM or JCMT) beam, to the observation s 
performed towards IRAS16293-2422 l|Parise et alJl2005h . 
As expected the simple conversion works reasonably well 
for the high lying lines, as they originate from the small in- 
ner envelope (small enough compared to the beam size to 
be approximated by a point). The accuracy of the con- 
version for the ground line at 464 GHz is less good as 
some of the emission originates from the outer envelope, 
which is affected by the convolution of the telescope beam. 
Nevertheless the simple conversion is still accurate to 15% 
as for this testcase, the HDO abundance is three orders of 
magnitude lower in the outer envelope compared to the in- 
ner envelope. The conversion would become critical when 
the abundance in the cold envelope is not low compared 
to the warm contribution. 

As a final remark, the plots of Fig. 01 [5] and H3 can 
be used with single dish observations without the need of 
convolution, for the signal, originating in the warm part of 
the envelope, is likely to be encompassed by the telescope 
beam. 



5.2. The low luminosity source case: L1448mm 

The previous predictions are of course dependent on the 
physical structure considered for the source. As an ex- 
ample, we also plotted the contours for a low-luminosity 
source, L 144 8mm, for which phys ica l characteristics ca n 
be found in l.Iorgensen et all l|2002|) or lMaret et al.ll|2004h . 
The results are presented in Appendix El The contours 
roughly look the same, except that the emission is weaker 
in L1448mm. 



5.3. Are space based telescope observations absolutely 
necessary towards low-mass protostars ? 

We saw in the previous sections that the observation of the 
studied lines from the ground allows to perfectly constrain 
the inner HDO abundance. However, the outer HDO 
abundance can only be constrained using lower energy 
transitions. Among the studied lines, only the 464 GHz 
can provide, in some specific conditions (namely when the 
outer abundance is not too low compared with the inner 
abundance), some information on x C oid- Unfortunately, the 
transitions that will be observable from space are at higher 
energy (see Fig.EJ. They will hence only provide informa- 
tion on the inner abundance, which is already correctly 
constrained with the transitions observed with ground 
based telescope. Moreover observations from the ground, 
when possible, have to be preferred owing to the better 
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Transition 


Frequency 


Telescope 


Conversion coetl 




IKAb 16293-2422 






{(jtiz ) 




1 K km s (1 m b Av J 


(Is. km s j 


(rv km s J 


(K km s j 








(erg s -1 cm -2 ) 


Conversion" 


Convolution 6 


Obs c 


ll,0 — 


80.578 


IRAM 


1.4xl0 -17 


0.19 


0.21 


0.40 


3l,2 — 22,1 


225.897 


IRAM 


4.1xl0~ 17 


1.7 


1.8 


1.7 


2i,i — 2i,2 


241.561 


IRAM 


4.4xl0~ 17 


1.5 


1.5 


2.0 


22,0 — 3l,3 


266.161 


IRAM 


4.9xl0~ 17 


0.20 


0.22 


0.21 


lo.i — Oo,o 


464.924 


JCMT 


4.0xl0~ 16 


6.3 


5.4 


5.5 



Table 2. Conversion factors of the line intensity in erg s 1 cm 2 into the velocity integrated main beam temperature 
T m hAv observed a t the IRAM and JCMT telescopes respectively, for the five lines obse rved towards IRAS 16293— 2422 
l|Parise et al.l2005j) . Last column quotes the T mb Av observed towards IRAS16293-2422 (|Parise et al.ll2005|) . a Predicted 
flux from the model simply converted into K km s _1 using the conversion factor. 6 Predicted flux derived using an 
accurate convolution of the signal with the telescope beam. 



spatial resolution. Table |3] lists the integration time re- 
quired to observe some selected transitions with a peak 
signal to noise of 5, at a resolution of 0.5km/s, with the 
various instruments IRAM 30 m, JCMT and HSO/HIFI. 
Note that in each case, only one receiver was assumed (sin- 
gle polarization mode) so that the integration time could 
in principle be lowered by a factor of \f2 if the observations 
are carried out in double polarization mode. An average of 
2 mm of water vapor has been assumed at the IRAM 30 m, 
where the average elevation of the source is 25° and 70° 
respectively for IRAS16293 and L1448mm. Weather band 
2 was assumed at JCMT, with a zenith angle of 40 and 
20 respectively. An efficiency factor of 0.125 was assumed 
for HIFI, as predicted in chopping mode. 

These predictions show that for sources as bright as 
IRAS 16293, observations can be carried out from the 
ground, and no Herschel time is thus necessary. But in 
the case of a low-luminosity protostar, where the same 
abundance profile of HDO was assumed, integration times 
start to become prohibitive from the ground. The 225 and 
241 GHz lines may remain observable, using double polar- 
isation receivers. Note that the integration time required 
with HIFI is prohibitive, the only reasonable transition 
being the 995 GHz. Thus, in the case of low-luminosity 
source, neither ground-based telescopes nor HSO/HIFI 
seem to be sensitive enough. The interferometer ALMA 
might provide in this case the unique possibility to con- 
strain the HDO abundance. 

6. Conclusions 

We reported theoretical predictions of the deuterated wa- 
ter line emission from the envelopes of low mass proto- 
stars. In this study, we have focussed only on the envelope 
emission, neglecting any HDO emission from the outflows, 
which is expected to be of minor importance based on 
theoretical and observational arguments. We have shown 
that the simultaneous observations of appropriately se- 
lected transitions permit to approximatively constrain the 
HDO abundance in the outer cold and inner warm parts 
of the envelope. The uncertainty on the evaporation tern- 



Transition IRAM 30m 


JCMT 


HSO/HIFI 


IRAS 16293 (30 L Q ) 


80 GHz 5.6 hours 






225 GHz 15 min 






241 GHz 30 min 






464 GHz 


29 min 




490 GHz 


21 min 


7 min 


509 GHz 


34 hours* 


104 min 


893 GHz 




16 min 


995 GHz 




5 min 



L1448mm (5 L ) 



80 GHz 


3500 hours 






225 GHz 


110 hours 






241 GHz 


227 hours 






464 GHz 




185 hours 




490 GHz 




132 hours 


73 hours 


509 GHz 




1562 hours* 


471 hours 


893 GHz 






59 hours 


995 GHz 






16 hours 



* The sensibility drops on the edge of the band. 

Table 3. Integration time required to detect the HDO 

lines with a peak signal over noise of 5 (see text). Note that 

only one receiver was assumed, so in principle those times 

can be lowered by a factor \/2 by using double polarization 

setup. 



perature mostly translates in an uncertainty on the inner 
abundance. The inner abundance can however be quite 
accurately determined using the 225 GHz line for which 
the emission has been shown to be independent from the 
evaporation temperature. 

The most important result of the present study is 
that, for bright low-mass protostars (L>10L Q ), observa- 
tions feasible from ground telescopes are enough to con- 
strain the HDO abundance profile and evaporation tem- 
perature, and no HSO observations are required. In the 
case of low- luminosity protostars (L<10L Q ), both present 
ground based telescopes and HSO seem to lack sensitivity 
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to observe HDO lines in a reasonable integration time. The van der Tak, F. F. S., Schilke, P., Miiller, H. S. P., et al. 2002, 

HDO abundance profile and evaporation temperature will A&A, 388, L53 

probably be directly constrained by future high resolution Vastel, C, Phillips, T. G., Ceccarelli, C, & Pearson, J. 2003, 

observations with the large sub/millimeter interferometer ApJ, 593, L97 

ALMA. 
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Appendix A: The predicted fluxes for the testcase 
of IRAS16293-2422. 

Table A.l. Predicted fluxes for the IRAS16293-2422 
testcase. 



Transition 


Freq. 


E up 




Line flux 




(GHz) 


(cm- 1 ) 


(s- 1 ) 


(erg s" cm" ) 



3.3,0 — 3 3 


1 


0.824671 


233.050 


5o,5 - 4 2 


2 


3.29676 


221.946 


4 3 ,i - 4 3 


2 


5.70278 


295.677 


22,0 — 2 2 


1 


10.2782 


109.269 


3 2 ,i - 4i 


4 


20.4600 


157.064 


53,2 — ^3 


3 


22.3077 


374.409 


32,1 — 32 


2 


50.2363 


157.064 


43,1 — 52 


4 


61.1859 


295.676 


6o,6 — 52 


3 


69.5506 


306.312 


li.o - ll 


1 


80.5783 


32.4939 


5i, 5 - 4 2 


2 


120.778 


225.862 


6i,e - 5 2 


3 


138.531 


308.612 


4 2 ,2 - 4 2 


3 


143.727 


221.832 


3 2 ,i - 4 


4 


207.111 


157.060 


3l,2 — 22 


1 


225.897 


116.456 


2i,i — 2i 


2 


241.562 


66.1781 


52,3 — 4 3 


2 


255.050 


303.989 


22,0 — 3i 


3 


266.161 


109.262 


52,3 — 52 


4 


310.533 


303.987 


62,5 — 53 


2 


313.751 


384.867 


3 3 ,i - 4 2 


2 


335.396 


233.016 


53,2 — 61 


5 


356.836 


374.402 


lo,i - Oo 





464.925 


15.4975 


33,0 — 42 


3 


479.947 


233.039 


3i,2 — 3i 


3 


481.780 


116.449 


2o,2 — li 


1 


490.597 


46.1612 


li,o — lo 


1 


509.292 


32.4844 


22,0 — 3o 


3 


537.793 


109.256 


2i,i — 2o 


2 


599.927 


66.1706 


3i,2 — 3o 


3 


753.411 


116.444 


4i,3 - 4i 


4 


797.487 


182.965 


4i, 3 - 3 2 


2 


827.263 


182.964 


2l,2 - li 


1 


848.962 


58.1067 


li,i - Oo 





893.639 


29.7880 


4 3 ,i - 5i 


4 


912.605 


295.656 


2o,2 — lo 


1 


919.311 


46.1517 


4i, 3 - 4 


4 


984.138 


182.961 


3o,3 — 2i 


2 


995.411 


91.3064 


2i,i - li 





1009.94 


66.1608 


52,3 — 5l 


4 


1161.95 


303.968 


4 2 ,2 - 4i 


3 


1164.77 


221.809 


5i,4 — 5i 


5 


1180.32 


265.208 


32,1 — 3i 


2 


1217.26 


157.036 


3i,3 — 2i 


2 


1267.04 


100.361 


2l,2 — lo 


1 


1277.68 


58.0972 


22,0 — 2i 


1 


1291.64 


109.239 


5i,4 — 5o 


5 


1297.81 


265.206 


3o,3 — 2o 


2 


1353.78 


91.2989 


32,2 ~ 22 


1 


1392.92 


155.357 


32,1 — 22 





1432.88 


157.032 


5i, 4 - 4 2 


3 


1444.83 


265.202 


4o,4 - 3i 


3 


1491.93 


150.121 


3 3 ,o - 4i 


3 


1500.99 


233.016 


3i,2 — 2i 


1 


1507.26 


116.426 


22,1 — 2i 


2 


1522.93 


108.890 


61,5 — 61 


6 


1615.63 


362.469 


3i,3 — 2o 


2 


1625.41 


100.353 


32,2 — 3l 


3 


1648.80 


155.350 



2.12x10" 


12 


2.6 


xlO" 


28 


4.34x10" 


13 


3.2x10" 


27 


4.17x10" 


10 


2.1x10" 


25 


3.63x10 


-9 


1.3x10" 


22 


9.58x10 


- 9 


2.6 


xlO" 


22 


1.64x10 


-8 


1.7x10" 


23 


2.08x10 


-7 


1.8 


xlO" 


20 


2.66x10 


-7 


1.3x10" 


21 


3.48x10 


- 9 


1.2x10" 


22 


1.32x10 


-6 


2.7x10" 


18 


1.33x10 


-6 


3.8 


xlO" 


19 


1.43x10 


-6 


9.7 


xlO" 


20 


2.80x10 


-6 


2.6 


xlO" 


19 


1.14x10 


-7 


4.2 


xlO" 


20 


1.31x10 




7.0 


xlO" 


17 


1.18x10 


-5 


6.5 


xlO" 


17 


1.78x10 


-5 


1.0 


xlO" 


18 


1.75x10 


-5 


1.0 


xlO" 


17 


1.78x10 




1.1 


xlO" 


18 


3.75x10 


-0 


1.1 


xlO" 


18 


2.61x10 




1.2 


xlO" 


18 


3.53x10 


-7 


5.8 


xlO" 


21 


1.69x10 


-4 


2.5 


xlO" 


15 


7.28x10 


-0 


4.6 


xlO" 


18 


4.74x10 




2.6 


xlO" 


16 


5.25x10 


-4 


6.4 


xlO" 


15 


2.32x10 


-3 


1.8 


xlO" 


15 


9.94x10 


-7 


1.8 


xlO" 


18 


3.45x10 


-3 


2.9x10" 


15 


5.90x10 


-3 


4.2x10" 


15 


1.31x10 


-4 


3.2x10" 


16 


1.28x10 


-3 


4.2x10" 


15 


9.27x10 


-4 


1.0x10" 


14 


8.35x10 


-3 


1.1x10" 


14 


2.74x10 


-6 


2.1x10" 


19 


1.56x10 


-3 


1.5x10" 


14 


1.07x10 


-2 


5.4x10" 


15 


7.04x10 


-3 


2.3x10" 


14 


1.56x10 


-3 


1.5x10" 


14 


2.19x10 


-2 


1.6x10" 


15 


2.00x10 


-2 


3.0x10" 


15 


2.97x10 


-4 


2.3x10" 


16 


1.89x10 


-2 


4.4x10" 


15 


3.91x10 


-3 


2.8 


xlO" 


14 


2.19x10 


-2 


3.1x10" 


14 


1.59x10 


-2 


6.3x10" 


15 


1.97x10 


-2 


5.4 


xlO" 


15 


5.30x10 


-3 


3.6 


xlO" 


14 


3.25x10 


-3 


7.3 


xlO" 


15 


3.54x10 


-3 


7.2 


xlO" 


15 


1.02x10 


-2 


8.6 


xlO" 


15 


3.07x10 


-2 


4.0 


xlO" 


14 


3.67x10 


-6 


7.9 


xlO" 


19 


6.58x10 


-3 


3.6 


xlO" 


14 


2.06x10 


-2 


1.1 


xlO" 


14 


5.84x10 


-4 


1.5 


xlO" 


16 


4.49x10 


-2 


5.4x10" 


14 


3.08x10 


-2 


1.1x10" 


14 
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continued from previous page 



Transition 


Freq. 

(GHz) 


E 

( c m 


- 1 ) 




A ui 




Line flux 
(erg s — cm 


4i 


4 


- 3i 


3 


1678 


58 


156 


343 


9 


92 x 10~ 


3 


4 


1x10" 


14 


61 


5 


- 60 


6 


1684 


Gl 


362 


468 


3 


47x 10~ 


2 


4 


3x10" 


15 


4o 


4 


- 3 


3 


1763 


56 


150 


115 


1 


20 x 10 


2 


5 


0x10" 


14 


4 2 


3 


- 4i 


4 


1818 


53 


217 


000 


4 


12 x 10~ 


2 


1 


0x10" 


14 


4 2 


3 


-3 2 


2 


1848 


31 


216 


999 


1 


06 x 10~ 


2 


1 


1x10" 


14 


4 3 


2 


-3 3 


1 


1872 


61 


295 


443 


6 


46 x 10~ 


3 


9 


6x10" 


16 


4 3 


1 


-3 3 





1877 


49 


295 


634 


6 


51 x 10~ 


3 


9 


7x10" 


16 


2-2 


1 


- 2o 


2 


1881 


29 


108 


883 


1 


21x10 


4 


7 


6x10" 


16 


3 2 


2 


- 3 


3 


1920 


43 


155 


344 


3 


04 x 10 


4 


1 


1x10" 


15 


4 2 


2 


- 3 2 


1 


1941 


80 


221 


791 


1 


23 x 10~ 


2 


1 


1x10" 


14 


5 3 


3 


-6! 


6 


1950 


15 


373 


620 


1 


1 1 x 10 




9 


4x10" 


19 


5 3 


3 


- 61 


6 


1950 


15 


373 


620 


1 


11 x 10~ 




6 


3x10" 


14 


5o 


5 


- 4i 


4 


1965 


55 


221 


900 


8 


28 x 10 


2 


4 


2x10" 


14 


4i 


3 


- 3i 


2 


1994 


29 


182 


937 


1 


66 x 10 


2 


4 


0x10" 


14 


4 2 


3 


- 4 


4 


2005 


18 


216 


995 


5 


46 x 10 


4 


8 


0x10" 


16 


5 3 


3 


- 60 


6 


2019 


13 


373 


618 


4 


52 x 10~ 


4 


3 


9x10" 


17 


5 2 


4 


- 5i 


5 


2031 


74 


293 


589 


5 


43 x 10 


2 


7 


9x10" 


15 


61 


5 


- 5 2 


4 


2064 


69 


362 


459 


4 


07 x 10 


2 


9 


5x10" 


15 


5i 


5 


-4i 


4 


2083 


04 


225 


816 


1 


98 x 10~ 


2 


3 


4x10" 


14 


4 3 


2 


- 5i 


5 


2087 


23 


295 


438 


9 


50 x 10 


6 


1 


6x10" 


IS 


5 3 


2 


- 5 2 


3 


2110 


99 


374 


360 


8 


12x10 


2 


3 


0x10" 


15 


5 2 


4 


- 5 


5 


2149 


22 


293 


587 


8 


51 x 10~ 


4 


4 


5x10" 


16 


5o 


5 


- 4 


4 


2152 


20 


221 


896 


2 


23 x 10 


2 


4 


0x10" 


14 


4 3 


2 


- 5 


5 


2204 


71 


295 


436 


3 


39 x 10 


4 


6 
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Appendix B: The low luminosity case : L1448mm 

We present here the results for the case of the low- 
luminosity source L1448mm. T he density and tem perature 
profiles where derived bv I J0rgensen et all (|2QQ2j) . 
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Fig. B.l. HDO line fluxes predicted for the case of the 
low-luminosity source L1448mm, for an assumed HDO 
abundance of 10" 10 and 10 -7 in the cold and warm parts 
of the envelope respectively. The upper panel presents the 
transitions observable from the ground, the middle panel 
the transitions observable with HSO/HIFI and the lower 
panel the transitions observable with SOFIA. The fre- 
quencies of the transitions that are going to be studied 
in more detail are indicated in GHz. 




Fig. B.2. Line intensity as a function of the HDO abun- 
dance in the outer and inner regions of the envelope, ai co id 
and x warm respectively, for the eight HDO transitions at 
464, 893, 509, 80, 241, 225, 919 and 995 GHz respectively 
towards the low-luminosity source L1448mm. Line fluxes 
are in ergs -1 cm~ 2 . In these calculations the mantle evap- 
oration temperature is 100 K. 



Fig. B.3. Line intensity as a function of the inner HDO 
abundance and the evaporation temperature, for the eight 
HDO transitions at 464, 893, 509, 80, 241, 225, 919 and 
995 GHz respectively towards the low-luminosity source 
L1448mm. Line fluxes are in ergs" 1 cm~ 2 . In these calcu- 
lations the outer HDO abundance is 10 -10 . 



